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ABSTRACT 

Theoretical studies of the physical processes guiding the formation and evolution of galaxies 
and galaxy clusters in the X-ray are mainly based on the results of numerical hydrodynam- 
ical N-body simulations, which in turn are often directly compared to X-ray observations. 
Although trivial in principle, these comparisons are not always simple. We demonstrate that 
the projected spectroscopic temperature of thermally complex clusters obtained from X-ray 
observations is always lower than the emission-weighed temperature, which is widely used in 
the analysis of numerical simulations. We show that this temperature bias is mainly related to 
the fact that the emission-weighted temperature does not reflect the actual spectral properties 
of the observed source. This has important implications for the study of thermal structures in 
clusters, especially when strong temperature gradients, like shock fronts, are present. Because 
of this bias, in real observations shock fronts appear much weaker than what is predicted 
by emission-weighted temperature maps, and may even not be detected. This may explain 
why, although numerical simulations predict that shock fronts are a quite common feature 
in clusters of galaxies, to date there are very few observations of objects in which they are 
clearly seen. To fix this problem we propose a new formula, the spectroscopic-like tempera- 
ture function, and show that, for temperature larger than 3 keV, it approximates the spectro- 
scopic temperature better than few per cent, making simulations more directly comparable to 
observations. 

Key words: Cosmology: numerical simulations - galaxies: clusters - X-rays: general - galax- 
ies - hydrodynamics - methods: numerical 



1 INTRODUCTION 

Clusters of galaxies have great cosmological relevance due to the 
privileged role they play in the hierarchical scenario of cosmic 
structure formation. In fact they represent the largest cosmic struc- 
tures that had time to undergo gravitational collapse and virialize. 
This characteristics makes galaxy clusters a powerful tool and a 
fundamental ingredient in many tests for the determination of the 
main cosmological parameters, like the matter density Q.,„ and the 
normalization and shape of the power spectrum of primordial fluc- 
tuations (represented by the spectrum amplitude Og and spectral 
index n, respectively). 

From a theoretical point of view, the relevant quantity to mea- 
sure is the cluster virial mass; this - however - can be directly 
determined only through gravitational tensing measurements, and 
even so many systematic effects make such measurement not an 
easy one. Alternatively, virial masses can be estimated indirectly 
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from related quantities, as d one with X-ray or millimetri c (e.g. the 
Sunyaev-Zel'dovich effect; ISunvaev & Zeldovichlll973l observa- 
tions, which - under appropriate assumptions - link the cluster mass 
to its temperature. Theoretical arguments suggest in fact the exis- 
tence of the so called scaling laws, tight correlations between mass 
and other global quantities, in this case the mean cluster tempera- 
ture (Kaiser 1986). Such correlations, confirmed an d calibrated by 
the results of hydrodynamical simulations (see, e.g., 'Navarro et aj 
1995; Evrard et al. 1996; Brvan & Norman 1998; Eke et al. 199^ 
allow for instance to link the observed cluster temperature func- 
tion to the theoretically estimated mass function. Another, more 
detailed, method to infer cluster masses is based on the solution 
of the equation of hydrostatic equilibrium, where however an as- 
sum ption on the spatial temperatu re distributi on is required (see , 
e.g., lEvrard et ail 1 19961; iTormen et alJll997t iThomas etalJll998t 
iRasia et aL2004. 



These arguments illustrate how important it is to have a 
reliable description of the thermal structure of galaxy clus- 
ters. This description needs also be spatially detailed, as re- 
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cent observational data with high spatial and spectral resolu- 
tion suggest that clusters are far from isothermal, and show in- 
stead a number of peculiar thermal features, like cold fronts 
(.Markevi tch & et al. 2000; Vikhlinin et al. 2001; Mazzo tta et al. 
l200lt iMa rkevitch et al. 2001; Mazzotta et al. 2003), c avities 
jMcNamara et alj I200Q; .Fabian et al. 2000; Blanton et a'i]|200lL 
McNamara et alj 1200 ll: M azzotta et al. 2002; Marke vitch et alj 



2001; Johnstone et al. 2002; Heinz et al. 2002; iYoung et alj|2002t 
ganders & Fabian 2002; Smith et al. 2002), blobs and filam ents 
iFabian et all200ll;lMazzotta et all2002l ; lMazzotta et alj2003h . 

A theoretical interpretation of these observations clearly re- 
quires state-of-the-art hydro-N-body simulations, which can be 
used to extract realistic temperature maps and/or profiles. The com- 
parison between real and simulated data, however, is complicated 
by different problems, produced both by projection effects and by 
instrumental artifacts, like instrument response, sky background 
and instrumental noise. Finally, a further complication can arise 
from a possible mismatch between the spectroscopic temperature 
Tspec estimated from X-ray observations and the temperatures usu- 
ally defined in numerical results. In fact, while Tspec is a mean pro- 
jected temperature obtained by fitting a single or multi-temperature 
thermal model to the observed photon spectrum, theoretical models 
fully exploit the three-dimensional thermal information carried by 
gas particles and so usually define physical temperatures. 

For the above reasons, possible biases can naturally arise when 
comparing the thermal structure of observed and simulated galaxy 
clusters. Using a set of 24 hydrodynamic-simulated clusters to ob- 
tain simulated images with quality similar to that e xpected in real 
Chandra observations. iMathiesen & Evrarjj200lh found that the 
spectroscopic temperature is lower than the mass-weighted temper- 
ature by roughly 10 to 20 per cent (see also Mathiesen et al. 1999). 
They claim that the origin of this bias is the excess of soft X-ray 
emission due to small clumps of cool gas that continuously merge 
into the intracluster medium. 

In order to have a more realistic comparison with the spec- 
troscopic fits, a different definition of temperature was then in- 
troduced to analyse the results of numerical simulations, namely 
the emission-weighted estimate both including and excluding the 
appropriate cooling function (see references in the next section). 
However, even in this case there is a tendency to infer significantly 
higher temperature values when the clu ster has a complex thermal 
structure (see, e.g.. lGardini et alJl2004) . In this context, the main 
goal of our paper is to propose and test a new formula for the theo- 
retically estimated temperature. Our formula should make simula- 
tions more directly and accurately comparable to observations. 

The plan of the paper is as follows. In Section 2 we introduce 
the problems originated by the projection of the cluster gas tem- 
perature along the line of sight, both from an observational and a 
theoretical point of view: the different definitions of projected tem- 
perature used in the literature are here introduced. In Section 3 we 
discuss in detail the effects of fitting realistic spectra - produced by 
multi-temperature thermal models - by single-temperature models, 
considering the specific case of Chandra observations. In Section 
4 we propose a new analytic formula capable of approximating the 
cluster spectroscopic temperatures measured from Chandra and 
XMM-NewUm observations with an accuracy always better than 
few percent. In Section 5 we test the performance of this new rela- 
tion by using the outputs of a high-resolution hydrodynamical sim- 
ulation, and compare them to those obtained adopting the widely- 
used emission-weighted estimate. Finally in Section 6 we discuss 
our results and summarize our conclusions. 



2 PROJECTING CLUSTER GAS TEMPERATURES 
ALONG THE LINE OF SIGHT: THE PROBLEM 

Although the atmosphere of galaxy clusters and groups is often far 
from isothermal, it is common practice to assume for it a single 
temperature, refereed as projected cluster temperature, in order to 
identify the thermal structure of the system along the line of sight. 
This is routinely done both by X-ray observers and by theoreticians, 
who provide projected temperature estimates of real and simulated 
clusters, respectively. In the following we review the techniques 
used to produce projected temperatures from simulated clusters and 
discuss why they may differ from the estimates obtained through 
the spectral analysis of observed clusters. 



2.1 Temperature from hydro-N-body simulations: mass- and 
emission-weighted temperatures 

Hydrodynamical N-body simulations provide information on the 
density and temperature of each gas element (be it particle or grid 
point) within the simulation box. Generally speaking, workers in 
the field of numerical simulations derive projected temperatures by 
simply calculating the mean weighted value of the gas temperature 
along the line of sight: 



^WdV 



(1) 



where T is the temperature of a gas element, dV is the volume along 
the line of sight and W is the weight function. 

In literature there are essentially two classes of tempera- 
ture projection: mass-weighted and emission-weighted. The first 
uses the mass of the gas element as weighting function (i.e. 
W = m; see. e.g.. [Kangrt aU1994 ; .Bartelmann & Steinmetz.l99S 



IMathiesen & Evrarj200ll) : 

_ ^niTdV 



(2) 



This temperature definition was first introduced because of its rele- 
vant physical meaning: the total thermal energy of the gas is simply 
E oc mrmw Despite of this, there has always been a general concern 
on the fact that such definition would give temperature estimates 
which differ significantly from what an X-ray observer would de- 
rive via spectral analysis. The main reason for this discrepancy is 
that the X-ray emissivity Eg is proportional to the square of the gas 
density (eg <>= rather than oc «), so it is expected that the spec- 
troscopic temperature, based on the number of emitted photons, 
is determined more by regions at higher density than by those at 
lower density. In order to make the temperature estimates of sim- 
ulated clusters more similar to the spectroscopic one, a new kind 
of temperature was introduced, named emission-weighted, where 
the weighting function is proportional to the emissivity of each gas 
element (W = e^): 



jA{T)n^TdV 
jA{T)n2dV 



(3) 



with A(T) the coo ling function and n the gas density (see, e.g., 
lNavarroetalll995h . 

Most of the works in the literature use the so called 
bolometric cooling func t ion A(r) = \n ZydE °c ^/T (see, e.g., 
Brva n & NormaiJ Il998l ; iFrenk & et all Il999t iMuanwong et alJ 
i200lft . implicitly assuming bremsstrahlung (free-free) emission, 
which is the dominant mechanism at temperatures larger than 3 
keV. 
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If not otherwise specified, all tiie emission-weighted tempera- 
ture estimates reported in this paper will refer to this definition. In 
other theoretical papers a much simpler definition is used, that as- 
sumes A(r) = 1 (see, e.g., Kang et al. 1994; Mathieseii^^Evrard 
l200lh . This definition is also known as emission-mea sure-weighted 
tempe rature. Only recently some authors (see, e.g.. lBorgani et al] 
|2004) tried to improve the correspondence between simulations 
and observations made with specific X-ray observatories, by adopt- 
ing the cooling function integrated on the specific telescope energy 
band: A(r) = J^™" EeJE, rather than the bolometric one. 

In the following we will show that none of these temperature 
definitions are accurate approximations of the observed spectro- 
scopic temperature. 



2.2 Spectroscopic projected temperature 

From the viewpoint of an X-ray observation, the cluster gas temper- 
ature is obtained by a fit of a thermal model to the observed spec- 
trum. Measuring a projected temperature is thus equivalent to find- 
ing a thermal model with temperature Tspec whose spectral proper- 
ties are as close as possible to the properties of the projected spec- 
trum. 

Now, from plasma physics we know that, if the emitting gas is 
a single temperature thermal plasma, its spectrum can be written as 
a linear combination of continuum and line emission processes: 



(.line 



(4) 



where E is the photon energy. For high gas temperatures (T > 
3 keV) and/or low metallicity (Z << I, with Z in units of solar 
metallicity) the continuum emission dominates over the line emis- 
sion. Let us assume, for now, that plasma metallicity is zero, so 
that we can set e^-'"' — 0. In this case the emission spectrum can be 
written as: 



„cont 
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(5) 



where n,, is the gas electron density. The function Gc{Z, T,E) is the 
total Gaunt factor, which is the sum of the Gaunt factors of three 
main different continuum emission processes, namely free-free 
(//, thermal bremsstrahlung), free-bound (fb), and two-photon 
(2y). 

Regardless of the functional form of the Gaunt factor, it is self- 
evident that, from a purely analytic point of view, the total spectrum 
induced by two thermally isolated plasmas with electron density ni 
and «2 and different temperature Ti and T2 can no longer be de- 
scribed by a single-temperature thermal model, with any tempera- 
ture T^. In fact, 

n2G,(Z,ri,£)-i-exp(-Jj-)-f 



nlGc{Z,T2,E)^exp{-i^_):^niGc{Z,T3,E)^^exp{- 



(6) 



unless Ti —72- This consideration is quite interesting, as it tells us 
that in principle the spectroscopic projected temperature is not at 
all a well defined quantity. In fact, as the spectrum of any single- 
temperature model cannot completely reproduce the spectral prop- 
erties of a multi-temperature source, the inferred spectroscopic 
temperature is a quantity that in principle may depend on the effi- 
ciency of the X-ray detector used for the observation and on the en- 
ergy band used to fit the spectrum. This is extremely important, as 
it tells us immediately that a lot of attention must be paid when we 
compare observational temperatures with the emission-weighted 
temperatures defined earlier in S I2.1I A proper comparison between 



simulations and observations requires the actual simulation of the 
spectral properties of the cluste rs via an X-ray obse rvatory simula- 
tor, like, for example, X-MAS jiGardini et ali2004 . 



3 SPECTROSCOPIC PROJECTED TEMPERATURE 
FROM Chandra OBSERVATIONS 

In the previous section we argued that from a purely analytic point 
of view the spectrum of a multi-temperature thermal model can- 
not be reproduced by any single temperature thermal model. In 
the real world things may be a bit different as, even assuming the 
most favorable conditions, observed spectra are affected by at least 
the following factors: i) convolution with the instrument response; 
ii) Poisson noise; iii) instrumental and cosmic backgrounds. These 
factors all conjure to distort and, at some level, confuse the ob- 
served spectra. Consequently it can happen that, under some cir- 
cumstances, observed spectra produced by multi-temperature ther- 
mal sources may be well fitted by single-temperature thermal mod- 
els which have little to do with the real temperature, but neverthe- 
less are statistically indistinguishable from it. In the present section 
we want to address exactly this issue. In order to do that, we gen- 
erated a number of multi-temperature spectra and performed the 
standard fitting procedure using a single-temperature model. 

To simulate the spectra we used the command FAKEIT in the 
utility XSPEC ^ Amaud.1 996: see , e.g., Xspec User's Guide version 
1 1 .2.x: IBorman &^mau3l200ll '). This command creates simu- 
lated data from the input spectral model by convolving it with the 
ancillary response files (ARF) and the redistribution matrix files 
(RMF), which fully define the response of the considered instru- 
ment, and by adding the noise appropriate to the specified inte- 
gration time. To begin with, in the following subsections we will 
concentrate on observations made with Chandra in ACIS-S con- 
figuration. In particular, all simulated spectra will be obtained by 
using the ARF and RMF relative to the aim point of the chip ACIS- 
S3. Later, we will discuss the differences between ACIS-S3 and 
the ACIS-I Chandra and MOS and pn XMM-Newton detectors. 
In order to make the simulated spectra more similar to the actual 
observations, we multiply each spectrum by an absorption model 
(WABS in XSPEC) that accounts for the galactic absorption. In this 
paragraph we fix the equivalent hydrogen column density of the 
absorption model to wh = 10^" cm^^. 

Our main interest is in the possible spectral distortion induced 
by instrument response, rather than in the spectral uncertainties 
connected with the data statistics. For this reason we will assume 
that the observed spectra have no background, and will rescale ex- 
posure times so that the total number counts per spectrum is very 
large (A' 350,000). 

All simulated spectra presented in this section are fitted by an 
absorbed single-temperature MEKAL model using the statistic, 
and leaving Z, and T as free parameters. Spectral fits are done 
in the 0.3 — 10 keV energy band. Regardless of the quality of the fit, 
from now on we will call spectroscopic temperature Tspec the best 
fit temperature value resulting from the above fitting procedure. 
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4 P. Mazzotta et al. 



Mekol models T1=1.0 (keV) T2=10.0 (keV) Mekal models Tl=2.5 (keV) T2=10.0 (keV) Mekol models T1=5.0 (keV) T2 = 10.0 (keV) 




0,5 1 2 5 10 0,5 1 2 5 10 0,5 1 2 5 10 



channel energy (keV) channel energy (keV) channel energy (keV) 

Figure 1. Upper panels: spectra of the input absorbed two-temperature thermal MEKAL model (solid line). The higher temperature model (dot-dashed line) 
is T2 = 10 keV in all panels. The lower temperature model (dashed line) is T\ = \ keV, T\ = 2.5 keV, and Ti =5 keV for the left, middle, and right panels, 
respectively. In all panels we use nn = 10^" cm^^ and Z = 0. Middle panels: simulated Chandra spectra corresponding to the input model in the upper panels 
and best fit model with one single thermal MEKAL model. Lower panels: residuals (in units of a) with error bars of size one. 




Figure 2. Left panel: reduced x~ corresponding to the best fit model of a single-temperature thermal model to the simulated spectra of a two-temperature 
thermal model in the case Z = 0. Each spectral fit has w 500 d.o.f. Solid and dashed histograms indicate the fits for which the null hypothesis probability 
is larger and smaller than 5 per cent, respectively. Right panel: percentile difference between the calculated {T) emission- weighted and spectroscopic (fspec) 
temperatures as a function of the temperatures of the lower and higher model components (T\ and T2, respectively). 



3.1 Fitting two-temperature thermal spectra with 

single-temperature models: the zero-metallicity case 

We will start by considering the simplest possible model, namely 
the one in which the gas metallicity is Z = 0. In this case the source 
emission is given by the sum of the continuum spectra (see Eq.|5} 
relative to each plasma component with different temperatures. For 
didactical purposes we first compare the simulated spectrum of a 



two-temperature thermal source with the best fit from a single- 
temperature model. We will assume that the two components of the 
thermal source have the same emission measure {koc J nldV). Let 
us consider tiu'ee examples, named spectrum 1, 2, and 3 for conve- 
nience. All spectra are characterized by the same higher tempera- 
ture component, fixed at 72 = 10 keV. The lower temperature com- 
ponent is instead = 1 keV, Tj = 2.5 keV, and Ti —5 keV for spec- 
trum 1, 2, and 3, respectively. In the left, middle, and right panels of 
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Figure 3. As Fig.0 but for metallicity Z = 1. 




Figure 4. As Fig.|2] but for metallicity Z = 1. 



Fig.Qwe report the two-temperature thermal input models, the sim- 
ulated source spectrum, the best fit of the single-temperature ther- 
mal model, and the residuals in units of the temperature dispersion 
a for spectrum 1, 2, and 3, respectively. The figure clearly shows 
that the source spectrum I (T] = I keV) cannot be fitted by a single- 
temperature thermal model. Conversely, spectrum 3 (Ti = 5 keV) 
is well fitted by a single temperature thermal model and actually is 
statistically indistinguishable from it. It is interesting to note that 
also spectra with low-temperature components Ti < 5 keV can ac- 
tually be fitted quite reasonably by a single-temperature model: in 
fact, although spectrum 2 (Ti = 2.5 keV) shows some small de- 
partures from its best fit single-temperature model (see the middle 



panels of Fig.Q), the x^/d.o.f=526/513 indicates that the fit is still 
statistically acceptable. 

Starting from the three examples above we can explore the 
intervals of and T2 required for the source spectra to be well 
fitted by a single-temperature model. To do so, we produced a large 
set of two-temperature source spectra with different Ti and T2, and 
estimated the spectroscopic temperature for each. 

The first interesting result we find is the following: con- 
sistently with what found by Mathie sen & Evrard 1 2Q01) , when- 
ever the lower temperature component is at J > 2 keV, then al- 
most any two-temperature spectrum can be well fitted by a single- 
temperature model. This is clearly shown on the left panel of Fig.|2| 
where we report the reduced relative to the best fit model of 
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Figure 5. Percentile difference between the calculated {T) emission- 
measure-weighted T (with weight W °= «j) and the spectroscopic (Tspec) 
temperatures as a function of the lower and higher model component tem- 
peratures (Ti and T2, respectively). 

a single-temperature thermal model to the simulated spectra of a 
two-temperature thermal model. The different line styles indicate 
if the fit is statistically acceptable or not. In particular, solid and 
dashed histograms refer to the fits for which the null hypothesis 
has a probability > 5 per cent (statistically acceptable fit) or < 5 
per cent (statistically unacceptable fit), respectively. We also no- 
tice that the reduced is always very close to unity, except in few 
cases where the lower temperature component is at T ~ 2 keV and 
the higher temperature component is at T > 8 keV. It is important 
to say that, for lower temperature components at T < 2 keV, our re- 
sults show that a single-temperature model is generally not a good 
fit, unless the difference — | is small (< 0.5 keV). In partic- 
ular, the lower the temperature of the low-temperature component, 
the smaller should be difference 1 72 — Ti | in order for its spectrum 
to be compatible with a single-temperature model. 

A second important result we find is that the spectroscopic 
temperature Tspec is always lower than the bolometric emission- 
weighted temperature, in agree ment with reports by some authors 
(see, e.g., Mathiesen & Evrardl EoOlL iGardini etai]l2004 . This is 
clearly shown on the right panel of Fig. |2| where we plot the 
percentile difference between the emission-weighted and spectro- 
scopic temperatures as a function of the lower and higher temper- 
ature components of the model: we notice that the larger | r2 — Ti | , 
the larger this difference, with values as large as 60 per cent for 
r2 < 12 keV that can be even larger for higher T2. 

3.2 Fitting two-temperature thermal spectra with 
single-temperature models: the case Z = 1 

We now study the effect of the line emission on the spectroscopic 
temperature estimates. To do that we repeat the analysis of S 13.21 
but now assuming Z = 1. 

As done in the previous subsection, we first show the spec- 
tra for the three example models, namely spectra 1, 2, and 3. On 
the left, middle, and right panels of Fig. |3| we report the two- 
temperature thermal input models, the simulated source spectrum, 
the best-fit single-temperature thermal model, and the residuals in 
terms of o for spectrum 1, 2, and 3, respectively. By comparing 
Fig. Q with Fig. |3| we notice that the presence of metals in the 
plasma induces the formation of lines that may not be reproduced 
by the best-fit single-temperature model, as clearly shown in the 



middle panels of Fig. |3| Although the continuum of the source 
spectrum is very well reproduced by the best fit single-temperature 
model (see. Fig.0, the line emission is under-predicted at 1 keV 
and over-predicted at 7 keV (see residual plot in Fig.|3}- This sug- 
gests that single-temperature models are not suited to fit multi- 
temperature spectra. Nevertheless, when the low-temperature com- 
ponent is high enough, then a single-temperature model can still fit 
very well a two-temperature source spectrum, and again it becomes 
statistically indistinguishable from it (right panels of Fig.|3j- This 
is because the higher the temperatures, the less important the line 
contribution to the final spectrum; in such cases the source with 
Z = 1 becomes more and more similar to that with Z = 0. 

As done in § 13.11 we produced a set of synthetic two- 
temperature source spectra with different T\ and T2, and estimated 
for each the spectroscopic temperature. In the left panel of Fig. |4] 
we report the reduced relative to the best fit single-temperature 
model. As expected, we find that in this case the number of source 
spectra that cannot be fitted by a single-temperature model is 
higher. Nevertheless, all source spectra whose lower temperature 
component is smaller than 3 keV continue to be statistically indis- 
tinguishable from a single-temperature model. 

Is is worth noticing that also in the case with Z = 1 we find that 
the spectroscopic temperature is always lower than the emission 
weighted one, as clearly shown in the right panel of Fig. |4| We 
also notice that the observed percentile difference in temperature is 
higher compared with the case with Z = 0, although only by a few 
per cent. 

We also tested the goodness-of-fit of a much simpler emission- 
measure-weighted temperature function, i.e. the one obtained by 
assuming A(r) = 1 (see ij |2.U . The percentile difference between 
the spectroscopic and emission-measure-weighted temperature cal- 
culated in this way is shown in Fig. |5] Paradoxically, we find 
that this simpler temperature definition fares better than the other 
one, even if it still overpredicts Tspec with differences that, for 
r2 < 12 keV, can be as large as 40 per cent. We must add that a 
different definition of Tgwi which uses the cooling function inte- 
grated in the telescope energy band (A(r) = /|™« e£rf£), gives 
results in-between the two discussed here. These results demon- 
strate that none of the emission-weighted temperature functions so 
far used in the literature actually provides a good approximation 
of the projected spectroscopic temperature obtained directly from 
observations. 

For completeness, we conclude this section by reporting the 
effects of projection on the hydrogen column density and metal- 
licity estimates. To do that we compare the fitted nn and Z of the 
single-temperature model with the input values. Percentile differ- 
ences are shown in the left and right panels of Fig.|6| respectively. 
Interestingly, we find that both the measured iif] and Z are different 
from the input values. This difference, however, is always within 
10 per cent^. 

4 APPROXIMATE FORMULA TO ESTIMATE 
Chandra AND XMM-A^eH'ton SPECTROSCOPIC 
TEMPERATURES 

In the previous section we have shown that, under some circum- 
stances, a multi-temperature source spectrum can be fitted by a 

^ Notice: the histogram in the right panel of Fig. |5| where the percentile 
variation of «/y is smaller than 10 per cent corresponds to regions where the 
fit is not statistically acceptable. 
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Figure 6. Percentile difference between the value obtained by fitting a single-temperature model to a two-temperature source spectra and the input value. Left 
and right panels refer to the equivalent hydrogen column density na (input value n^/.^ = 10^" cm^^) and to the metallicity Z in solar units (input value Zi„ = 1). 




Figure 7. Percentile difference between the calculated (T) and the spectroscopic (rspec) temperature a a function of the lower and higher model component 
temperatures (Ti and T2, respectively). To calculate T we use the spectroscopic-like temperature Ts\ defined in EQ. I14l with the value of a chosen in order to 
minimize the value of A. Left panel: observational conditions defined as in S I3.1l (to be compared with the right panel of Fig.|2). Right panel: observational 
conditions defined as in S I3.2l (to be compared with the right panel of Fig. 0. 



single-temperature model, providing the estimate Tipec- We have 
also demonstrated that T^pec is always lower than the bolometric 
emission- weighted temperature Tew In this section we wish to de- 
rive a projected emission-weighted temperature formula that can 
better approximate the spectroscopic one. 

The idea is quite simple: given a multi-temperature thermal 
emission we want to identify the one temperature whose spectrum 
is closest to the observed spectrum. From now on, we will call 
this "spectroscopic-like" temperature. If we assume two thermal 
components with constant densities n\, 112, and temperatures Ti, 
T2, respectively, requiring matching spectra means that 



"K(2,7'l)^exp(~jf.)+n2^(Z,r2,)^exp(- 
«.A^(Z,7;,i)A_exp(-j|), 



E 



(7) 



where A is an arbitrary normalization constant and ^(Z, T) is a 
parametrization function that accounts for the total Gaunt factor 
and partly for the line emission. 

Both Chandra and the XMM-Newton are most sensitive to 



the soft region of the X-ray spectrum, so we can expand both sides 
of Eq.0in Taylor series, to the first order in E/kT: 



«?C(z,ri)^(i-J;)+.2^(z,r2)^(i-J,; 
«.AC(z,7;,)^(i- 



E 

kTj, 



(8) 



By equating the zero-th and first-order terms in E, we finally 
find the equation defining the spectroscopic-like temperature: 



«2^(z,r,)/r;/V«2c(z,r2)/r2' 



1/2 



(9) 



n\^{Z,Ti)lTll^+nl(,{Z,T2)lTll^' 

The extension of Eq.|9|to a continuum distribution of plasma 
temperatures in a volume V is trivial: 



^n\{ZJ)lT^I'^dV 
jn^(,{ZJ)/T^I-dV' 



(10) 



To calculate the function ^(Z, T) we consider the fact that it 
gives, together with w^, the relative contribution of the spectral 
component with temperature T to the total spectrum. The ampli- 
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Table 1. Value of a and corresponding minimum A for different observa- 
tional conditions 
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Detector 
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a 
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0.46 


1.5 
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MOS 
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2.3 


5 
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34 


6 


MOS 


1 
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0.91 


4.0 
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MOS 
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1.0 


0.55 


2.2 


8 


MOS 


1 


1 


1.0 


0.96 


24 



tude of this normalization contribution is set partly by the total 
Gaunt factor, partly by the line emission. In principle, their rela- 
tive importance depends on the details of the instrument/analysis 
used (e.g. the energy band used for the fit) and on the observational 
conditions (e.g. low or high galactic absorption regions). However, 
for temperatures and metallicities realistic in clusters of galaxies 
(0 < Z < 1, 0.5 keV < r < 20 keV) the total Gaunt factor de- 
pends essentially on the temperature, and can be approximated by 
a power-law relation: 

G,(Z,^,£)«.(^/keV)^ (11) 

where r| ~ 0.25 (see, e.g.. lMewe & GronenschilJl98lh . The effect 
of the lines is to increase the plasma emissivity, especially at low 
energies: this effect can also be approximated by a power-law of 
the temperature. Thus in the following we will assume for ^(Z, T) 
the following functional form: 

^(z,r)oc(r/keV)", (12) 

where the parameter a depends on the specific observational con- 
ditions and on the used instrument, so it needs to be appropriately 
calibrated. We do this by adopting the following procedure: we se- 
lect some specific observation conditions, simulate a set of two- 
temperature source spectra as explained in § |3| and calculate the 
mean value of the percentile variation of T^i with respect to Tspec: 




In the previous formula the sum is extended to all the A' simulated 
source spectra with 3 keV < Ti < 20 keV and 3 keV <T2 <20 keV. 
The variable a is obtained though a minimization procedure of A. 

To show how well Tg\ reproduces the actual Tspec. we apply 
it to the cases discussed earlier in 13. H and i; 13.21 The minimiza- 
tion procedure of A returns a = 0.86 and a = 0.65 for the Z=0 
and Z=l, respectively. The resulting percentile variation between 
T^l and Tjpec for these two cases is shown in Fig.Q on the left and 
right panels, respectively. The figure clearly shows that, unlike Tew, 
Tgi provides a good estimate of the observed spectroscopic temper- 
ature, at a level better than 2-3 per cent. 

We now study how the power index a depends on the observa- 
tional conditions for all the Chandra and XMM-Newton detectors. 
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Figure 8. Value of a and coiTesponding minimum value of A for the ob- 
servational conditions reported in TablelTI Open circles and filled triangles 
refer to Chandra ACIS-S3 and XMM-Newton MOS, respectively. The er- 
ror bars show the range in a where A < 5 per cent (see text). The horizontal 
line indicates a = 0.75. 



and how it changes as a function of the minimum energy con- 
sidered in the spectral fitting. In TableQwe report for some of the 
cases considered the value of a obtained using the minimization 
procedure, and the corresponding value of minimum A. The cases 
reported in the table are only a few examples of those we actually 
considered. They represent extreme situations for galaxy clusters, 
in the sense that we explicitly considered strong and null galac- 
tic absorption («// = 10^^ cm^^ and n// = cm^^, respectively), 
and solar (Z=l) and zero (Z=0) metal abundances. The values of a 
and the corresponding minimum A from TableQare also shown in 
Fig.m Open circles and filled triangles refer to Chandra ACIS-S3 
and XMM-Newton MOS, respectively. From Fig.|8|it is clear that, 
as expected, the value of a minimizing the discrepancy between 
Tsi and Tspec depends on the actual observation conditions. Never- 
theless, we notice that: i) a and A are very close to each other for 
Chandra and XMM-A'ewton; ii) the minimum A is generally very 
low and smaller than 4 per cent. Following these considerations we 
explored the interval in a for which the corresponding value of A 
is smaller than 5 per cent. The results are reported in Fig.|8| where 
this range is shown as y-axis error bars. We notice that all the error 
bars overlap. This is important as it means that, regardless of the 
observational conditions, we can select a value of a (which we call 
a,,) so that can reproduce Tspec with an accuracy better than 5 
per cent on average. As shown by the solid horizontal line in Fig.|8| 
a good choice for Chandra and XMM-A'ewton is a,, = 0.75. 

The accuracy of Ts\ when we adopt a„ = 0.75 is shown in 
Fig. |9| and Fig. 1101 where we report the percentile variation be- 
tween the values of and Tspec obtained in different observa- 
tional conditions with the Chandra ACIS-S3 detector. In Fig. |9| 
and Fig. llOl the values of Tspec were obtained by fitting the spectra 
down to iinjin = 0.3 keV and = 1 keV, respectively. The four 
panels in each figure correspond to («// = x lO^^cm^^, Z = 0), 
{nn = 1 X lO^^cm"^^ 2 = 0), (nn = x lO^^cm"^, Z = 1), and 
(nn = 1 X lO^^cm^^, Z = 1). As expected, for most of the tem- 
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Figure 9. Percentile difference between tlie calculated {T) and tlie spectroscopic (Tspcc) temperatures as a function of tlie lower and higher model component 
temperatures (T^ and T^, respectively). To calculate T we use the spectroscopic-like temperature defined in Eg. 1141 assuming a = a„ = 0.75. The panels 
from up to down, left to right, conespond to the observational conditions 1 to 4, reported in TablelTlwith ACIS-S3, respectively. 



perature combinations of the source spectra, T^\{'a.u = 0.75) gives 
a good estimate of Tspec at a level better than 5 per cent, while the 
difference is higher only for a few spectra (see, e.g., upper-right and 
lower-left panels of Fig.|9|and Fig. llOt . However, the only spectra 
for which the discrepancy is higher than 5 per cent are those on the 
border of the temperature plane, i.e. those whose the lower tem- 
perature component in the source model is T = 3 keV. As already 
explained in S 13. H and S 13.21 this temperature represents a sort of 
"borderline" value for Chandra and 'KMM.-NewUm. Thus, the rea- 
son of such higher discrepancy is simply explained by the inability 
to properly define a spectroscopic projected temperature for these 
spectra, due to the poor quality of the fit with a single temperature 
model (see left panels of Fig.|2|and|4}. 

To summarize the results of this section, we claim that the 
spectroscopic-like temperature defined by: 



Til (a) i 



/n2r«/r3/2rfV/' 



(14) 



with a = a„ = 0.75, gives a good approximation of the spec- 
troscopic temperature Tspec obtained from data analysis of Chan- 
dra and XMM-A'ewron observations. 

We notice that Eq.^|can be rewritten in the general form of 
Eq.0 



Tsl - 



w 



fWdV 
n- 

7-3/4 • 



(15) 



It is interesting to note that weights each thermal component di- 
rectly by the emission measure but, unlike Tgw, inversely by their 
temperature to the power of 3/4. This means that, beside being bi- 
ased toward the densest regions of the clusters, the observed spec- 
troscopic temperature is also biased toward the coolest regions. In 
the next sections we will discuss some important implications of 
this bias. 



5 TESTING SPECTROSCOPIC-LIKE TEMPERATURE 
ON HYDRO-N-BODY SIMULATION 

Using a simplified two-temperature thermal model we have shown 
that the spectroscopic-like temperature provides a good approx- 
imation to the spectroscopically derived temperature Tspec- Here we 
wish to extend this study by showing the goodness of in repro- 
ducing Tspec in a more general case. 

In order to do that, we follow the same procedure used in 
iGardini et alJ i2004h . where we use the output of an hydro-N-body 
simulation and the X-ray MAp Simulator (X-MAS) to compare 
Tspec with Tew for a simulated cluster of galaxies. Here we use the 
same simulation to extend the comparison to T^\. 

Let us briefly remind that the cluster simulation used in 
lGardinieta'i]<200 4 was selected from a sample of 17 objects ob- 
tained by re-simulating at higher resolution a patch of a pre-existing 
cosmological simulation. The assumed cosmological framework 
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Figure 10. As Fig.|9] but the panels from up to down, left to right, correspond now to the observational conditions 5 to 8 reported in TablefTlwith ACIS-S3, 
respectively. 



is a cold dark matter model in a flat universe, with present mat- 
ter density parameter = 0.3 and a contribution to the density 
due to the cosmological constant flyv = 0.7; the baryon content 
was set to iig = 0.03; the value of the Hubble constant (in units 
of 100 km/s/Mpc) is /; = 0.7, and the power spectrum normal- 
ization is given by Og = 0.9. The re-simulation method we used, 
called ZIC (for Zoomed I nitial Conditions), is described in de- 
tail in iTormen et alJ il997 l). while an extended discussion of the 
properties of the whole sample of these simulated clust ers is pre- 
sented elsewhere iTormen et alj|2004 liasia et alJl2004) . Here we 
remind only some of the characteristics of the cluster used in this 
paper. The simulatio n was obtained by us ing the publicly avail- 
able code GADGET jSoringel et ai]|200ll) ; during the run, start- 
ing at redshift zj^ = 35, we took 51 snapshots equally spaced in 
log(l + z), from z = 10 to z = 0. The cluster virial mass at z = 
is 1.46 X 10'^/?^'Mq, corresponding to a virial radius of 2.'ihr^ 
Mpc; the mass resolution is 4.5 x X^h^^Mr.. per dark particles and 
5 X per gas particles; the total number of particles found 

inside the virial radius is 566,374, 48 per cent of which are gas 
particles. The gravitational softening is given by a 5hr^ kpc cubic 
spline smoothing. 

In the left panel of Fig. II II we report the emission- weighted 
temperature map of the simulated clusters overlaid to the clus- 
ter flux distribution from Gardini et al. (2004). The figure clearly 
shows that the cluster is far from isothermal. Among other features, 
we point to the reader the presence in this map of two shock fronts. 
The first is in the lower-left corner, and is produced by the motion 



toward the cluster centre of the innermost of the two subclumps 
present in that region. This front has a post-shock gas temperature 
of approximately 1 8 keV. A similar, but weaker, shock front is lo- 
cated in front of the subclump, at the centre of the Eastern side of 
the cluster. For comparison, on the right panel of Fig. ll ll we report 
the spectroscopic-like temperature map of the same cluster, over- 
laid to the cluster flux distribution. An immediate thing to notice is 
that the map of appears cooler than the map of Tew This is con- 
sistent with the fact that, unlike Tew, Ta is biased toward the lowest 
values of the dominant temperature components along the line of 
sight. Another very important point is that both shock fronts, which 
are clearly evident in the emission-weighted temperature map, are 
no longer detected in the map. This aspect will be further dis- 
cussed in §|5|below. 

In the following subsections we will compare Tspec, Tew and 
in two cases of practical interest: the determinations of the pro- 
jected cluster temperature map and the determination of the radial 
profile. To do that we will use the 300 ks^ Chandra ACIS-S 3 X- 
MAS "observation" of this cluster from iGardini et all f2004h . As 
explained in that paper, the observation for this particular cluster 
was performed by fixing the cluster metallicity to Z = 0.3Zp; and 
the column density to Nu = 5 x 10^^ cm^^; data were analysed by 
applying the stan dard procedur es and tools used for real observa- 
tions. We refer to lGardini et alj f2004) for more details on the data 
analysis. Here we just remind that T^pec is obtained by fitting the 
spectra in the [0.6,9.0] keV energy band with a single/temperature 
absorbed MEKAL model, after fixing the cluster redshift, metallic- 
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Figure 11. Comparison of the map of the emission-weighted temperature Tew (left-panel) and the spectroscopic-like temperature Ts\ (right-panel) for the 
simulated cluster of galaxy. Both maps are obtained by using the gas particles of the hydro-N-body simulation and are binned in 1" pixels. The contour levels 
con'espond to the cluster flux distribution. The circle shows the cluster region from which the temperature profile sispayed in Fig. llSl has been extracted. 



ity, and hydrogen column density to the values used as inputs to 
compute the Chandra observation. 



5.1 Cluster projected temperature map 

As mention ed above, we used th e spectroscopic temperature map 
published in lGardinietal1j2004 . That map was obtained by sub- 
dividing the cluster image in squares large enough to contain at 
least 250 net photons. The projected spectroscopic temperature 
map Tspec is shown in Fie. 1121 In order to make a direct comparison 
of this map to those with emission-weighted and the spectroscopic- 
like temperatures, we decreased the resolution of the latter to match 
the resolution of the former. 

In the left panel of Fig. 1131 we report the same map of Tew 
shown on the left panel of Fig. 1111 but re-binned as the map of 
Fig.EH 

To highlight the temperature differences between the spectro- 
scopic and the emission-weighted temperature maps, in the right 
panel of Fig. 1131 we show the percentile difference of (Tew — 
?spec)/7spec- For better visualization we only show the pixels where 
this difference is significant to at least 3o confidence level, i.e. 
\{Tew/~ 7spec)/Ospec| > 3, being Ospec the 68 per cent confidence 
level error associated with the spectroscopic temperature measure- 
ment. This plot clearly shows that there are many regions where 
the difference between Tew and Tspec is significant to better than 
> 3o confidence level. Furthermore, for these pixels the discrep- 
ancy ranges from 50 per cent to 200 per cent and even more. Of 
particular relevance are two cluster regions showing a shock front 
in the emission-weighted map, in the left panel of Fig.^3 The map 
of the percentile difference between Tgw and Tspec shows discrep- 
ancies of 100-200 per cent, indicating that shock fronts predicted in 
the emission weighted map are no longer detected in the observed 
spectroscopic map. 

Conversely, in the left panel of Fie. ll4l we show the same map 
of Ts\ shown in the right panel of Fie. 1111 but re-binned as the spec- 
troscopic temperature map of Fie. 1121 As we did for the emission- 
weighted temperature map, in the right panel of Fie. 1141 we show 
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Figure 12. Spectroscopic temperature map of the simulated cluster of 
galaxy as derived from the spectroscopic anal ysis of the Ch andra "ob- 
servation" with the package X-MAS (from Gardini et all2004l) . 



the percentile difference of {T^\ — T^^^) /T^^ec- Again we only visu- 
alize the pixels for which the difference is significant to at least 3cj 
confidence level. The presence in this map of fewer pixels clearly 
indicates that the match between and Tspec is much better than 
the one between Tspec and Tew. Furthermore, most of these pix- 
els shows very small temperature discrepancies, being smaller than 
10 per cent. Only in 6 pixels we find a slightly higher temperature 
discrepancy, but in any case smaller than 20 per cent. 

This demonstrates that the spectroscopic-like temperature 
gives a much better estimate of the observed spectroscopic temper- 
ature than the widely used emission-weighted one. It is worth not- 
ing that in the previous map 3 out of 6 pixels where the discrepancy 



12 P. Mazzotta et al. 




Q s t5 ;p .5C i.no i^o jqo 

T,i«ev) (■[„ -T^J/r^ (S3 



Figure 13. Left panel: emission-weighted temperature map of the simulated cluster of galaxies shown in the left panel of Fig. II li re-binned to match the 
spatial resolution of the spectroscopic temperature map shown in Fig. 1121 Right panel: percentile difference between the spectroscopic and emission-weighted 
temperature maps. In this map we show only the regions where the significance level of the temperature discrepancy is at least 3c, i.e. | (Tiw — Tspe^ ) / Ospcc I > 3, 
where Ospcc is the 68 per cent confidence level error associated to Tspec ■ 




Figure 14. Left panel: spectroscopic-like temperature map for the simulated cluster of galaxies shown in the right panel of Fig. ll ll re-binned to match the spatial 
resolution of the spectroscopic temperature map shown in Fig. II 21 Right panel: percentile difference between the spectroscopic and the spectroscopic-like 
temperature maps. In this map we show only the regions where the significance level of the temperature discrepancy is at least 3a, i.e. | {T^\ — T,pec)/cJspec| > 3, 
where Ospec is the 68 per cent confidence level error associated to r^pec . 



is between 10 and 20 per cent correspond to cluster regions with 
very low surface brightness. We believe that the observed discrep- 
ancy in this case is simply related to the very poor statistics used to 
determine Tipec- In the other 3 cases the relatively higher discrep- 
ancy is instead related to the fact that the thermal components of 
these regions that contribute to the spectra have a quite large spread 
of temperatures and, because the lower (dominant) thermal compo- 
nent is at r < 3 keV, Tspec cannot be unequivocally identified: a 
proper spectral analysis would require a fit with a two-temperature 
model. 



It is very important to say that, unlike the emission-weighted, 
the map on the right panel of Fig. ^| does not show big discrep- 
ancies between and Tspec in both shock cluster regions. This 
clearly indicate that does a much better job than Tgw in pre- 
dicting the projected spectral properties of such peculiar thermal 
features. 
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Figure 15. Temperature profiles of the simulated cluster of galaxies. Filled 
circles refer to the mean spectroscopic temperature profile from Gardini 
et al. (2004) and correspond to the region shown as a circle in Fig. 1111 
error bars are at 68 per cent confidence level for one interesting parame- 
ter. Solid and dotted histograms refer to the mean emission-weighted and 
spectroscopic-like temperature profiles directly extracted from the simula- 
tion. 

5.2 Cluster projected temperature profile 

To conclude this section we also tested the accuracy of in 
predicting the observed spectroscopic temperature profile. We 
again used the spe ctroscopic projected temperature profile of 
lGardinietal1<2004 . This was obtained by extracting spectra from 
circular annuli centred on the cluster centre, out to the radius iden- 
tified by the circle in Fig. 1111 The size of the bin was chosen 
in order to have approximately the same number of photons in- 
side each aimulus. The spectroscopic temperature profiles Tspec, to- 
gether with their relative 68 per cent confidence level errors Ospeci 
are shown as filled circles in Fig. 1151 In the same figure we show the 
emission-weighted temperatu re and the spectrosco pic-like temper- 
ature. As already discussed in iGardini et al]<2004 . we notice that 
the emission-weighted temperature Tew profile does not reproduce 
the spectroscopic temperature profiles Tspec- In particular we con- 
firm that the spectral temperatures are systematically lower than the 
emission-weighted ones. Conversely, the spectroscopic-like tem- 
perature profile provides a much more accurate estimate, and falls 
within the error bars of the "observed" spectroscopic temperature 
profile everywhere but in two annuli. 



6 DISCUSSION AND CONCLUSIONS 

In this paper we have studied the problem of performing a proper 
comparison between temperatures obtained from the data analysis 
of X-ray observations (i.e. projected spectroscopic temperatures, 
Tspec) and temperatures derived directly from hydro-N-body sim- 
ulations (projected emission-weighted temperatures Tgw)- In 12.21 
we show analytically that Tspec is not a well defined quantity. In 



fact, it results from the fit of a single-temperature thermal model 
to a multi-temperature source spectrum; however, since the former 
cannot accurately reproduce the spectral properties of the latter, it 
follows that Tspec cannot be unequivocally identified. Generally 
speaking, this means that a reliable comparison between simula- 
tions and observations can only be done through the actual sim- 
ulation of the spectral properties of the simulated clusters: these 
properties will then be directly compared with the observed ones. 
Nevertheless, observed spectra are affected by a number of factors 
that distort and confuse some of their properties. In some circum- 
stances multi-temperature thermal source spectra may appear sta- 
tistically indistinguishable from a single-temperature model. In §|3| 
we study this aspect focusing our attention on observations of clus- 
ters of galaxies made using the CCD detectors of Chandra and 
XMM-A'ewton. These detectors are characterized by having both 
a similar intermediate energy resolution and a similar energy re- 
sponse. From our study we find two very important results that for 
convenience we summarize below. 

(i) Given a multi-temperature source spectrum, if the lowest 
dominant temperature component has Ti > 2 — 3 keV, then a fit 
made with a single-temperature thermal model is statistically ac- 
ceptable regardless of the actual spread in temperature distribution. 

On the other hand, multi-temperature sources with Ti <2 keV 
will most likely require multi-temperature spectral models. This is 
equivalent to say that Tspec can be properly defined only for spec- 
tra with Ti > 2 — 3 keV. For lower temperatures the identification 
of multi-temperature observed spectra with a single temperature is 
not appropriate. It is important to say that this result is intrinsi- 
cally related to the characteristics of the X-ray detector used for 
the observation (e.g. its energy pass band and resolution) and not 
to a possible inadequate photon statistics. In fact our result was ob- 
tained in the limit of high spectral photon number (see §|3j. It is 
self-evident that this conclusion will be true also for all the detec- 
tors whose spectral properties are similar or worse than those on 
board of Chandra and XMM-Newton, while this will not be true 
for X-ray spectrographs with either a much larger energy pass band 
or a much higher energy resolution. 

(ii) The emission-weighted temperature Tgw, originally intro- 
duced to provide a better comparison between simulations and 
observations, in practice does not properly estimate the spectro- 
scopic temperature Tspec- In particular Tew tends to overestimate 
Tspec- This mismatch depends on the thermal inhomogeneity of 
the observed multi-temperature source: the larger the spread in 
temperature of the dominant components in the observed spec- 
trum, the larger the discrepancy. In §|4]we derive a new formula, 
the spectroscopic-like temperature function T^i (see Eq. I15> . and 
show that Tsi can approximate Tspec to a level better than 10 per 
cent regardless of the temperature spread. It is worth noticing that 

weights each thermal component directly by the emission mea- 
sure but inversely by their temperature to the power of 3/4. This 
explains why the observed Tspec is biased toward the lower values 
of the dominant thermal components. 

These two results have important observational and theoreti- 
cal implications for the study of X-ray clusters of galaxies. First 
of all the equivalence between observed "hot" multi-temperature 
spectra and single-temperature models implies that, using the CCD 
detectors of Chandra and XMM-A'ewton, it is observationally im- 
possible to disentangle the single or multi-temperature nature of 
any observed spectrum whose projected temperature is higher than 
3 keV by performing a simple overall X-ray spectral analysis. In 
addition, the temperature derived from this spectral analysis is not 
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Figure 16. Tempera ture profiles of the galaxy cluster 2A 0335 (from 
iMazzotta et aljl20O3l) . Filled triangles and open circles refer to the pro- 
jected cluster temperature profile in the northern and southern sectors, re- 
spectively. Open squares indicate the overall projected cluster temperature 
profile obtained from a circular radial analysis. The upper and the lower his- 
tograms show the profiles for the emission-weighted and spectroscopic-like 
temperatures, respectively. These histograms have been obtained by com- 
bining the northern and southern temperature profiles. 



the emission-weighted value, but it is biased toward the lowest 
dominant thermal component of the overall spectrum. The con- 
sequences for studies of temperature profiles in clusters are im- 
mediate. In virtually all works on the subject, in fact, the cluster 
temperature profile is derived by extracting spectra from concen- 
tric circular or elliptical annuli centred on the cluster X-ray peak. 
Together with the cluster gas distribution, this temperature profile 
is used to estimate the cluster mass by assuming hydrostatic equi- 
librium. It is self-evident that, if the cluster gas temperature distri- 
bution is azimuthally asymmetric, the temperature profile derived 
from the radial analysis is biased toward lower temperature val- 
ues and so does the estimated mass. To better show this aspect we 
discuss the results of the data analysis of the cluster of galaxies 
2A 0335-1-096 selected, as example, among the many published in 
the literature. In particular, a ll the temperature mea surements dis- 
cussed below are taken from^ Mazzotta et aP J2003h . In Fig. 1 161 we 
compare the projected temperature profiles of the cluster of galax- 
ies 2A 0335-1-096 extracted from different sectors. Filled triangles 
and open circles indicate the temperature profiles obtained from the 
Northern (from -90° to 90°; angles are measured from North toward 
East) and Southern (from 90° to 270°) sectors of the clusters. As 
highlighted by Mazzotta et al. ( 2003), the temperature profile of 2A 
0335-1-096 is clearly azimuthally asymmetric. We notice that in the 
100-200 arcsec radial interval the temperature profile of the south- 
ern sector is between 30 and 50 per cent hotter than the Northern 
one. The open squares in Fig. 1161 indicate the cluster temperature 
profile obtained from the circular radial analysis. For convenience 
we report as histogram the expected temperature profile obtained 
by combining the Northern and Southern temperature profiles us- 
ing the emission-weighted formula. From this figure it is evident 
that, consistently with what discussed so far, the circular radial tem- 



perature profile is significantly lower that the expected emission- 
weighted profile and is biased toward the values of the Northern 
profile. For completeness in the same figure we added as histogram 
the temperature profile obtained by combining the Northern and 
Southern temperature measurements using the spectroscopic-like 
formula. This latter profile is perfectly consistent with the mea- 
sured one proving once more that, unlike the emission-weighted 
temperature, our spectroscopic-like formula provides a very good 
approximation to the spectroscopic temperature measurements. 

It is important to say that, by construction, the difference be- 
tween Tew and mainly depends on the complexity of the pro- 
jected thermal structure of the cluster gas. It is clear that the find- 
ings of this study are especially relevant for major-merger clusters 
rather than for relaxed clusters. They may have also important im- 
plications for the study of all structures with very strong tempera- 
ture gradients like, for example, the shock fronts. In fact, because of 
the temperature bias of Tspec, shock fronts in real observations ap- 
pear much weaker than the predictions of virtually all the emission- 
weighted temperature map published in literature. This has been 
shown in 15. H and is evident in Fig. II II From this figure we imme- 
diately see that the two shock fronts clearly visible in the emission- 
weighted map are no longer detected in the spectroscopic-like tem- 
perature map. This temperature bias may explain why, although 
simulations predict that shock fronts are quite common in clusters 
of galaxies, to date we have very few observations of clusters in 
which they are clearly present. 

To conclude we stress once more that the emission-weighted 
temperature function may give a misleading view of the actual gas 
temperature structure as obtained from X-ray observations. Thus, 
since the emission-weighted temperature has no physical meaning 
(unlike the mass-weighted one), here we propose to theoreticians 
and N-body simulators to finally discard its use. We remind that 
great attention must be paid when comparing simulations and X- 
ray observations. Under the most generic conditions such compar- 
ison can only be done through the actual simulation of the spectral 
properties of the simulated clusters, so that software packages like 
X-MAS I Gardin i et al. 2004) become fundamental. Nevertheless, if 
the cluster temperature is sufficiently high and if the spectral prop- 
erties of the detector used for the observation are similar or worse 
that the ones of Chandra and XMM-Newton we showed that our 
proposed spectroscopic-like temperature function may be consid- 
ered an appropriate tool for the job. 
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